Early stages of the shear banding instability in wormlike micelles 
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We study the early stages of the shear banding instability in semidilute wormlike micelles using 
the non-local Johnson-Segalman model with a two-fluid coupling of the concentration <j> to the shear 
rate 7 and micellar strain W_. We calculate the "spinodal" limit of stability for sweeps along the 
homogeneous intrinsic flow curve. For startup "quenches" into the unstable region, the instability in 
general occurs before the homogeneous startup flow can attain the intrinsic flow curve. We predict 
the selected time and length scales at which inhomogeneity first emerges. In the "infinite drag" 
limit, fluctuations in the mechanical variables 7 and W_ are independent of those in <j>, and are 
unstable when the slope of the intrinsic flow curve is negative; but no length scale is selected. For 
finite drag, the mechanical instability is enhanced by coupling to 4> and a length scale is selected, 
in qualitative agreement with recent experiments. For systems far from an underlying zero-shear 
demixing instability this enhancement is slight, while close to demixing the instability sets in at low 
shear rates and is essentially demixing triggered by flow. 

PACS numbers: 47.50.+d Non-Newtonian fluid flows- 47.20.-k Hydrodynamic stability— 36.20.-r Macro- 
molecules and polymer molecules 



Control of morphology and stability is vital for pro- 
cessing many complex fluids, e.g. polymeric, liquid crys- 
talline, and surfactant fluids. While much is understood 
close (or relaxing) to equilibrium, strongly driven sys- 
tems suffer in comparison. For many complex fluids, the 
intrinsic constitutive curve of shear stress E versus shear 
rate 7 is non-monotonic. For semi-dilute wormlike mi- 
celles, theory predicts the form ACEG of Fig. [j] Q, f§. 
In the regime of decreasing stress, 7 c i < 7 < 7 C 2, homo- 
geneous flow is unstable || Q and the system splits into 
bands of different shear rates 7 = 7^,7h jj R pi U, M, 
with a steady state flow curve ABFG. Constitutive mod- 
els augmented with interfacial gradient terms have cap- 
tured this behaviour j^, |l2). In this paper we study 
the linear stability of initially homogeneous states. We 
also predict the length and time scales at which inho- 
mogeneity emerges after a startup "quench" , i. e. impos- 
ing a shear rate in the unstable regime. Our analysis is 
the counterpart, for this driven phase transition, of the 
Cahn-Hilliard (CH) calculation for demixing in undriven 
systems. 

Shear startup experiments j|, [H| reveal (i) a 
metastable regime 7^ < 7 ~ 7 c i of slow approach to 
the banded state, and (ii) an unstable regime, onset at 
7 ~ 7 c i, where the stress can massively overshoot £ se i 
before subsiding rapidly to it. In Ref. |Q, the overshoot 
coincided with concentration fluctuations that emerged 
perpendicular to the shear compression axis at a selected 
length scale 0(l/im). The authors attributed this to 
the Helfand-Fredrickson (HF) pi coupling of flow to 



concentration. This features in two component systems 
with well separated relaxation times (e.g. polymer and 
solvent) jl§ [L7| |T|, |||, §(J. The slow component 
(polymer) tends to migrate to regions of high stress. 
If the plateau modulus increases with polymer concen- 
tration, positive feedback enhances concentration fluc- 
tuations and can shift the spinodal of any nearby CH 
instability @. 




FIG. 1: Schematic flow curve. 



Here we model the initial stage of the unstable kinetics 
using the non-local Johnson-Segalman (d-JS) model |2^| 
within a 2-fluid framework ^) that incorporates con- 
centration dynamics. The 2-fluid model considers sepa- 
rate force balance equations for the micelles (velocity v m , 
volume fraction <fi) and solvent (velocity v_ s ). These are 
added to give overall force balance for the average veloc- 
ity v = <t>v m + (1 — <ft)v s (Eq. |l|) and subtracted for the 
relative velocity v rcl = v m — v s , which in turn specifies 
the concentration fluctuations (Eq. 0): 



2 



P (9t + u.V) v = D t v = V.G(0) K ~ <t> V 
</> 2 (l-0) 2 



2 V.(l>T ]m D m + 2V. (1 - 0) ?? S; D a - Vp, 



A0 = -Y-0(i-</>h 



ei = -V- 



Y-G((t>)K ^5F 2V- 4>VmR° m 2V- (l-0)?7 s £ o s 



(1) 
(2) 



G{4>)W_ is the viscoelastic micellar backbone stress, W Q /3 = -g# • -g^ — <5 Q/ 3 is the local strain that would have to be 
reversed to relax this stress, and G is the plateau modulus. The free energy F comprises osmotic and elastic parts: 
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where /" is the osmotic susceptibility and f the equilibrium correlation length for concentration fluctuations. The 
Newtonian stress 20?y m i2^ n describes fast micellar processes (e.g. Rouse modes) with Lr the traceless symmetric 
micellar strain rate tensor, and 2(1 — (t>)VsD^ s is the solvent stress. Incompressibility determines the pressure p. 
The drag coefficient £ (Eq. ||) ensures that the force Cn rc i impedes relative motion. The micellar diffusion coefficient 
D oc f"/C. We have omitted negligible inertial corrections to Eqs. ([!]) and (||) |pTf . 

Eq. ([l]) is the Navier Stokes equation generalised to include osmotic stresses. Eq. (||) is a generalised CH equation 
in which micelles diffuse in response to gradients in the osmotic force V^F 1 /^] an d m the viscoelastic stress; for 
dG/d4> > (assumed here), HF feedback occurs. The micellar strain obeys d-JS dynamics 63]: 



(dt+v m - Y)W = a{D -W + W-D ) + (W-Q -ft -W) + 2D - 



W 



I 2 



T(<t>) T{<f>) 



V 2 W 



(4) 



where 2ft_ m = Vu m -(Vu m ) T with (Vv m ) Q/ 3 = 9 a (wm)/3- 
t(4>) is the Maxwell time; I is a length that could, for 
example, be set by the mesh-size. The slip parameter a 
measures the fractional stretch of the micelles compared 
to the flow. For \a\ < 1 (slip) the intrinsic flow curve 
^xyij) is capable of the non-monotonicity of Fig. [|. 

We study planar shear between infinite plates at y = 
{0,L} with (w,Vf),VA») in the (i, y,z_) directions. At 
the plates we assume d y 4> = d y 4> = 0, d y W_ = 0, and 
no slip. For controlled strain rate (assumed through- 
out) Jq dyj(y) — const. Unless stated, we use model 
parameter values at <fi — 0.11 from rheological data 
for CTAB(0.3M)/NaNO 3 /H 2 O ||, and light scattering 
(DLS) data for CTAB/KBr/H 2 (§|; we calculate the 
drag as £ = 67r^ -2 [g3| where fj = 4>n m + (1 — 4>)r) s - We 
extrapolate G((j)), r((f>) } Z?(0), C(^)> to (j) < 0.11 using 
scaling laws for semidilute wormlike micelles. We fix a 
(assumed independent of (j)) by comparing with Cates' 
model 0; and use units in which G(4> = 0.11) = 1, 
t(<P = 0.11) = 1 and L = 1. The homogeneous intrinsic 
flow curves T, xy (j, (j)) that satisfy d t v = d t <f) = d t W_ = 
are shown as dashed lines in Fig. |2|. The region of 
negative slope ends at a "critical" point (f> c w 0.015: 
CPCl/NaSal in brine || shows the same trend. 

We linearise in fluctuations J2 k [Sj, 8W, S(j)]^ e l -- +ut 
about these homogeneous states. The stability analysis 
ujk,aY.k,a = M. k Y.k,a gives the normal modes v k a (where 
a is the mode index), each encoding a set of relative 
amplitudes of £7, 5Wij , 5<f>. The lower spinodal lies where 



the largest branch oj k of the dispersion relation first goes 
positive (unstable) as the background homogeneous state 
is swept up the intrinsic flow curve from 7 = 0; the upper 
spinodal is defined likewise, for sweeps down from large 
7. We consider only k = ky; the stability of k — kz is 
unaffected by shear in our model. 

In the limit of infinite drag C — * 00 at fixed V ■ | , 
the coupling of £7 and S W to 8(f) is disabled. Fluctuations 
in the mechanical subspace [£7, SW xy , SW XX , SW yy ] then 
obey conventional (uniform-</>) d-JS dynamics, in which 
any homogeneous shear state in the regime of negative 
constitutive slope d'Sxy/dj < is unstable: the spinodal 
is given by circles in Fig. ||. The concentration fluctu- 
ations independently obey conventional (zero-shear) CH 
dynamics, with a demixing instability when D < 0. We 
consider only flow-induced instabilities here, D > 0. 

For finite drag, HF feedback couples these subspaces 
and can induce instability even if dTixy/dj > and D > 
0. For model parameters from the data of Refs. |m| p5| , 
the mechanical instability is enhanced slightly by con- 
centration coupling (squares in Fig. [^). This effect in- 
creases near a zero-shear demixing instability: see the 
diamonds and triangles in Fig.|^, for smaller D((f> = 0.11). 
(The lobe of instability at high shear rate is discussed in 
Ref. plj .) Within some simplifying assumptions |5l[ , a 
qualitative expression for the lower spinodal is D d * y + 
§ W^f% = ^ where 22 = (a - l)W xx + (a + l)W yy < 0, 
C<xC and £ = f 

Enhancement of flow instabilities by positive feedback 
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FIG. 2: Dashed lines: intrinsic flow curves for <j> = 0.11, 0.091, 
0.072,0.053,0.034,0.015 (downwards). Symbols: spinodals for 
the uncoupled limit f — ► oo (o); coupled model with D(<f> — 
0.11) taken from DLS (□); coupled model with artificially 
reduced D (0,A). Inset: zoom on large 7. 



with concentration was first predicted by the remarkable 
insight of Schmitt et aZ.j2(|. However they directly as- 
sumed a chemical potential /i = ^(7). Although this 
is equivalent to our approach in the limit of adiabatic 
stress response (assumed in [^6|), below we find that the 
dynamics inside the spinodal are dictated by the micel- 
lar stress response. The spinodal is unaffected since the 
dynamics are adiabatic here by definition: the above con- 
dition corresponds to Schmitt's Eq. (24). 

Experimentally the spinodal is found via sweeps along 
the intrinsic flow curve, though in practice banding can 
occur prematurely via metastable kinetics ]l3| ] . The same 
ambiguity arises in the spinodal of conventional fluid- 
vapour demixing, defined via quasistatic compression. 

We now study the early-time kinetics in the unsta- 
ble regime. We cannot assume that the system starts 
in a homogeneous state on the unstable intrinsic flow 
curve, since the Maxwell time t needed to prepare such 
a state is longer than the typical time scales l/u>k of 
the instability itself (except very near the spinodal) [p8| , 
and so most startup flows go unstable before the intrinsic 




FIG. 3: Instability occurring as the homogeneous startup flow 
evolves towards the intrinsic unstable constitutive curve. 

flow curve can be attained. When the rheometer plate 
is first set moving the shear rate rapidly homogenises, 



on the Reynolds timescale pL 2 /i] <C r. The shear stress 
E(t) = GW(t) + 777 then starts to evolve towards the 
intrinsic flow curve, with Wit) initially given by the ho- 
mogeneous startup solution of Eq. (|J). Without noise, 
it would attain this flow curve at some time r ss = 0(t). 
But of course the system has noise. The stability problem 
is now time-dependent: L0k,a(t)vf. tOl (t) = M, (t)v fc a (t) 
since M =M k (7, Wit) ,4>). At first all dispersion branches 
are stable. Then at some time to < r ss a branch 
goes unstable. The size of the growing fluctuations at 



t > t is A k (t) 



exp 



rion for detectability is logA — O(10). 



A rough crite- 
which defines a 



k— dependent time scale Ti ns %(k) via J^""*^ dt'u>k(t') — 
O(10). In most regimes, fluctuations emerge fastest at a 
selected wavevector k* , because of a peak in the disper- 
sion relation Wfe(t) (Fig. ||). We thus define the overall 
time scale of the instability to be r inst = r inst (fc*) (Fig.||). 

Fig. [| gives the unstable dispersion branch Wfc(t) and 
the eigenvector at the maximum in u>k for three startup 
flows. Ti ns t is marked in each case. (We used a continu- 
ous k axis though in practice only harmonics k = nir are 
allo wed.) For the pure mechanical instability (£ — > 00, 
Fig. 4(a)), Uk is cut off at large k by interfaces and at 
low k (beyond the gap size) by Reynolds effects. Between 
there is a plateau where uik is set by the viscoelastic s tress 
response, with no clearly selected lengthscale. Fig. [i(b)| 
is for a coupled model far from a zero-shear CH demixing 
instability (with the spinodal given by circles in Fig. ||). 
Concentration coupling has enhanced the mechanical in- 
stability at short length scales; in competition with the 
interfacial terms this selects a length scale as seen 

experimentally [pi. The eigenvector at k* (Fig. 4(b), 
bottom) is mainly still in the £7 and S W directions, as ex- 
pected. The plateau of the mechanical instability is still 
apparent at smaller k, since this shear rate would have 
been unstable even for £ — > 00. Concentration diffusion 
cannot keep pace at these length scales, and is absent 
from the eigenvector. Fig. 4(c) is for a system close to 
zero-shear demixing, for which the lower spinodal is at a 
shear rate far below the spinodal of the pure mechanical 
instability (triangles in Fig. |2|). Our imposed shear rate 
is just inside this spinodal, so the mechanical plateau is 
absent from the dispersion relation, leaving only the dif- 
fusive concentration-coupled branch. The eigenvector is 
dominated by 6cj>, so this instability is essentially demix- 
ing, triggered by flow. It occurs slowly enough that the 
intrinsic flow curve is reached before fl uctu ations grow 



appreciably: Tj nst is off the scale of Fig. 4(c). 

In summary, we have studied the early time kinet- 
ics of the shear banding instability in the d-JS model 
with 2-fluid coupling to concentration. The spinodal on- 
set of pure mechanical instability with negative consti- 
tutive slope dTixy/dj < (mathematically indistinguish- 
able from a flow-induced transition in, for example, ne- 
matic liquid crystals ^7|) is shifted by coupling to con- 
centration. For startup quenches deep in the unstable 
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FIG. 4: Startup dispersion relation (top) and eigenvector at max. in Uk (bottom). Z = \{a - l)W xx + |(1 + a)W y 
Y= 1(1-0)^-1(1 + 0)^. 



region the instability in general occurs before the homo- 
geneous startup flow can attain the intrinsic flow curve. 
An initial length scale is selected only if the instability is 
coupled to concentration. These results are qualitatively 
consistent with recent experiments on wormlike micellar 
solutions |2l. 
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